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Abstract We calculate magnetic field fluctuations above a conductor with a nonlocal response (spatial dispersion) and 
consider a large range of distances. The cross-over from ballistic to diffusive charge transport leads to reduced noise 
spectrum at distances below the electronic mean free path, as compared to a local description. We also find that the mean 
free path provides a lower limit to the correlation (coherence) length of the near field fluctuations. The short-distance 
behavior is common to a wide range of materials, covering also semiconductors and superconductors. Our discussion 
is aimed at atom chip experiments where spin-flip transitions give access to material properties with mesoscopic spatial 
resolution. The results also hint at fundamental limits to the coherent operation of miniaturized atom traps and matter 
wave interferometers. 

PACS. 42.50.Ar Photon statistics and coherence theory; 42.50.Nn Quantum optical phenomena in absorbing, am- 
plifying, dispersive and conducting media; cooperative phenomena in quantum optical systems; 72.10.-d Theory of 
electronic transport; scattering mechanisms; 74.25. N- Response to electromagnetic fields; 42.50.Lc Quantum fluctua- 
tions, quantum noise, and quantum jumps 
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Electromagnetic fluctuations have been playing a key role in 
physics ever since Planck discovered the black-body spectrum. 
They have universal properties at distances from a body large 
compared to the thermal (Wien) wavelength. In fact, the noise 
spectrum is telling a lot about material properties in the near 
field, due to the links provided by the Kirchhoff law and the 
fluctuation-dissipation theorem 1 1 1. In practical applications like 
magnetic resonance imaging, this near-field noise is a limiting 
factor for detecting biological signals, for example |2,3 1. In the 
case of a metallic body, electromagnetic fluctuations depend 
mainly on the conductivity and can reveal details about charge 
transport in the bulk. We focus here on a linear cuiTent-field 
relation in the general nonlocal form (spatial dispersion) 



i{r,u!) = X! / d^^'<^mn{r,r',uj)En{r',uj) 



(1) 



where the dependence on both r and r' contains the crossover 
from ballistic to diffusive in the motion of charge carriers. This 
introduces the mean free path iasa characteristic length scale. 
The conductivity (or dielectric) tensor now depends on both 
frequency and wave-vector in Fourier space (spatial disper- 
sion) pHTO) . Further nonlocal effects are introduced by the de- 
tails of the surface and the surface scattering of charge carri- 
ers |^T,Ta-i6|. 

The anomalous skin effect is a famous consequence of the 
nonlocal bulk response. It is typically discussed in the regime 
of high frequencies, where the classical skin-depth 5{uj) falls 



below the mean free path and does no longer describe screening 
of magnetic fields correctly ||4|[17][18) . Other relevant physical 
phenomena are Thomas-Fermi (Debye-Hiickel) screening and 
Landau damping |19J, connected to plasma screening due to 
mobile charges in the metal and electron-hole generation (in- 
ternal photo effect), respectively. One can generally expect a 
reduction of field fluctuations close to the surface of a nonlo- 
cal metal as compared to local theory, because the bulk fields 
are better screened and escape less easily into the surround- 
ing space. The nonlocal response of surfaces has also been dis- 
cussed for other observables outside the range of the anomalous 
skin effect. Its implications have been worked out for the dis- 
persion relation of surface plasmon modes p3p6) , the transfer 
of heat via near field radiation f20"21'|, the van der Waals inter- 
action across a electrolyte |22, 23], and the Casimir interaction 
between two metallic half-spaces p4] - p6) . 

A motivation for this paper is the observation that spin-flip 
transitions of ultracold atoms held in miniaturized chip-based 
magnetic traps (atom chips) are sensitive to magnetic fluctua- 
tions in the near field of a metal L27j-j30J. Here, atoms are prob- 
ing surface properties at somewhat exotic frequencies in the 
radio or microwave band, much below the visible to ultravio- 
let frequency range of conventional metal spectroscopy. (Note, 
however, that the microwave band is routinely used in super- 
conductor experiments.) At the same time, the coiTesponding 
wavelengths are in the micron range because atoms illuminate 
the surface with their near field. Modes characterized by these 
frequencies and wave vectors (parallel to the surface) lie in 
the evanescent sector, way below the light cone; in particular 
their in-plane wavevector p is not restricted as in propagat- 
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ing vacuum fields. Therefore values p ^ Xji cannot be ex- 
cluded where spatial dispersion is clearly relevant. Previous 
work 1 31-341 has considered spin-flip transitions for the mi- 
crotrap scenario in the local limit, identifying them as relevant 
challenges to miniaturization below the micron scale. We also 
mention the results of Ref. |35 1 on nuclear spin relaxation that 
are reproduced and generalized here, including the skin effect 
and covering a wider range of distances. 

Another motivation is the study of spatial correlations of 
thermal near field radiation. These differ strongly from the black- 
body limit where the wavelength provides a universal correla- 
tion (or coherence) length (36) . For an overview on the coher- 
ence of thermal radiation see Ref. | |37j . Electric field correla- 
tions were discussed previously in Refs. pS ■ 40 1 for homoge- 
neous media and in Refs. |4T-'46l for the near field of bod- 
ies. Surface charge and current correlations have been stud- 
ied in the high-temperature limit in Refs. | p3l[47|[48| . In the 
electric case, the mean free path t did not emerge as a charac- 
teristic length scale, neither in the distance dependence of the 
noise spectrum nor in the spatial autocorrelation function [44 1. 
This may be related to sum rules and efficient screening at the 
surface. The magnetic field behaves differently because at dis- 
tances around the mean free path there is a crossover in the 
noise spectrum |35|. We show here that the field correlations in 
a plane parallel to the metal surface become more coherent at 
short distances (z <C (), and that the correlation length involves 
the mean free path. 

This work is organized as follows. In Sectionl2]we outline 
the calculation of magnetic spectra and spin-flip rates and give 
an overview on the length scales and effects that have an im- 
pact on these quantities in a system with nonlocality. Then, a 
specific nonlocal model for the bulk and surface response is in- 
troduced and used to obtain the near field asymptotes. We also 
briefly consider materials other than metals. In Sec.|3]field cor- 
relation functions above local and nonlocal metals are analyzed 
with respect to their correlation length. Sec. HI summarizes the 
main results and reflects on their relevance for experimental 
setups. Two Appendices give details on the electromagnetic 
Green's tensor near a surface and the calculation of the non- 
local reflection coefficients. 



2 Magnetic noise and spin-flip losses 

2.1 Noise spectra and spin flips 

Our main quantity of interest is SfAv^uj), the spectral den- 
sity per unit frequency of the magnetic field cross-correlation 
(i, i — X, y, z). This spectrum can be calculated with Green's 
function techniques, as outlined in Appendix [A] According to 
the fluctuation-dissipation theorem. 



S^At, uj) = 2hn{uj) Im U^jir, r, uj) 



(2) 



where n{ijj) is the Bose-Einstein distribution and T-Lij the mag- 
netic Green's tensor defined in Eq.( 30 1. 



Recall that the Green's tensor gives the field radiated by a 
pointlike dipole source. In the presence of a surface, it there- 
fore splits in two terms, the first one being the same as in free 
space [Eq. ([3T]l], the second one describing the magnetic field 



reflected from the surface. In the situations considered here, the 
latter term dominates the spectrum f49l. For example, at a fre- 
quency of 100 MHz and at room temperature, the spectrum of 
a surface at a distance of 1 /xm exceeds the black body spec- 
trum by 15 orders of magnitude. The free space term in % can 
therefore be safely neglected and it is sufficient to consider the 
reflected Green's tensor, which is conveniently expressed in the 
Weyl representation as a two-dimensional Fourier integral 



'Hij{z,u}) = ^ / dppn 
^7^ Jo 



rsi^,p) + -^^-^rp{uj,p) ] X 

C K 



x[5ij - ZiZj] + 2—rs{uj,p)ziZj 



(3) 



Here, s and p label the two principal polarizations, n ~ p — 
ixP' jt? is the propagation constant in vacuum, and z the unit 
normal to the surface. Details on the reflection coefficients rg, 
Tp are given in AppendixJA] 

Most of this work willconsider near field noise, where large 
values of the perpendicular wave vector p 3> i^/c (evanescent 
waves) dominate the response and nonlocal effects become rel- 
evant. In this regime, the p-polarization involving r^ is sup- 
pressed by the prefactor us^ J^ck)^ in the integrand of Eq. (|3]l. 
The analysis can thus be restricted to s -polarization for our pur- 
poses (magnetic field vector in the plane of incidence). 

The magnetic noise spectrum has been measured via the 
loss rate of atoms from modern chip-based atom traps 1 27 ■ 30) . 
An expression for the atomic transition rate due to fluctuations 
of the magnetic field can be obtained from Fermi's Golden 
Rule p8) or a master equation approach 1 32 34) 



..ah ..ha 



^^-w = E^^I^^.^(r 



-'^ah 



)• 



(4) 



Here a (b) labels a magnetic sublevel that is trapped (not trapped) 
in the static magnetic field of the atom chip, nf^ ~ (a|/ii|&) is 
the matrix element of the magnetic dipole operator, and cuab 
the resonant Bohr frequency. For magnetic moments in the or- 
der of a Bohr magneton fis, the prefactor in Eq. (|4]) translates a 
spectrum of 130 pT^/Hz to a transition rate of one per second. 
Rates in this low range have been measured near conducting 
surfaces using ultracold atoms as a probe | 33|50pT) . The mag- 
netic near-field noise flips the spin of trapped atoms, leading to 
loss from a magnetic trap and setting a fundamental limit to 
the coherence in these setups p7l|28) . Conversely, this process 
offers a way of probing material properties. 

For atoms trapped in their electronic ground state, the mag- 
netic moment is dominated by the contribution of the electron 
spin. We evaluate the matrix elements in Eq. (pi for simplicity 
by ignoring the quantum numbers of the nuclear spin. We can 
then consider a two-level system of which one state is magnet- 
ically trapped. For a static trapping field in the zz-plane that is 
tilted by an angle 9 relative to the surface normal z, the mag- 
netic dipole matrix elements read 1 32 52) 



M- = ^(cos^, 



-1, sine 



(5) 



Above a planar surface, the noise correlations are diagonal [see 
Eq.d6]l below] so that the spin-flip rate is proportional to 2 + 
sin^. The more general case including hyperfine structure can 
be found in Ref. 
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the Green's tensor can be approximated by the scaling laws 



b) 



distance z [^m] 



Figure 1. a) Spin flip (loss) rate near a conducting half-space de- 
scribed by the nonlocal Boltzmann-Mermin conductivity (ISjfTTl. The 
local description [Drude conductivity l |12| )] and the asymptotic expres- 
sions of Eq. ^ (dashed) are shown for comparison. The length scales 
A — vf /i-Op, I = uf/t, and 5 [Eq. jTl] illustrate the Thomas-Fermi 
screening length, the mean free path and the skin depth. The parame- 
ters are for gold at T = 300K(cr = 2.5 X W Q.'^ra-^,l^ 18 nm, 
7 = 6.7 X 10^"^ s~^) and the surface impedance is calculated with a 
specular boundary condition [Eqs. l |35[|36^ ]. Spin flips are driven by 
fields at the Larmor frequency uj/2-k — 1 MHz and oriented parallel 
to the surface [6 — in Eq. jsl]. Losses due to the free space black 
body spectrum are much smaller and not visible on this scale, 
b) Loss rates near gold surfaces with different purities. We vary the 
ratio 'j/iOp between relaxation rate and plasma frequency in the con- 
ductivity. The lowest curve coincides with Fig.flji). Note how the in- 
termediate regime i <^ z <^ S opens up in the dirty limit. The leftmost 
arrow marks the Thomas-Fermi screening length A — vp/ujp, while 
S is the normal skin depth HI. 



2.2 Overview: near-field noise 



Nonlocal effects can be expected to become visible on a length 
scale in the order of the mean free path £ of ballistic transport 
of charge carriers, as was akeady conjectured by Rytov and 
coworkers 153). Numerical calculations of the spin-flip rates 
for neutral atoms near a metal surface with and without a non- 
local response are shown in Fig.fTk) and b). Clearly, there are 
three different asymptotic regimes of the distance between the 
atom and the surface, two of which involve distances much 
larger than £, where the surface spectrum cannot be distin- 
guished from a local one. We shall find that in these regimes. 
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Here, the skin depth of the normal skin effect is given by 



6iL0) 



6{uj) < z 

£ < z < 5{uj) (6) 

0.077 1 , z < £ . 



(7) 



JIqCTUJ 



where a is the local limit of the DC conductivity [Eq. ( 12 1 be- 
low]. The local regime [first two lines in Eq.(|6]l] were given 
already in Refs. 1 31 32) . 

To make this qualitative behavior understandable, we pro- 
pose an interpretation in terms of an active surface volume: 

i) When z ^ <5(a;), the normal skin effect screens noise 
from deep in the bulk so that only a skin layer of thickness 
6{uj) contributes to the noise. The noise is proportional to the 
squared non-retarded fields ^ 1/r^ of current loops , integrated 
over the surface - this explains the power law l/z"* and the 
proportionality to the skin depth S{uj) in Eq. (J6|, upper line. 

ii) At smaller distances i -^ z <^ <^{^)^ ^ medium-filled 
half-sphere of radius ^ z effectively contributes to the noise. In 
addition, the probe particle now resolves individual current ele- 
ments rather than loops. The noise then arises from the squared 
fields of these current elements (^ 1/?"^), integrated over the 



volume of the half-sphere, as explained in Ref. 1 54 1 



The previous cases i) and ii) have been observed experi- 
mentally in the kHz to MHz range with sensitive magnetome- 
ters 1 551 and with trapped ultracold atoms |j50|. In this paper, 
we address the regime 

iii) of the extreme near field in a nonlocal conductor, z ^ 
£: the ballistic (rather than diffuse) motion of charge carriers 
creates spatial correlations in the fluctuating current field. This 
reduces the number of mutually uncorrected volume elements 
in the half-sphere introduced in ii) above, and hence lowers the 
noise power Note that the limiting value in Eq. (J6]) scales with 
l/((5^£) which is actually independent of the relaxation time 
in a Drude conductor The magnetic noise is related to Landau 
damping, or equivalently to the thermal excitation of electron- 
hole pairs flOj. This regime is therefore quite universal, and 
we show in Sees. 12.51 and 12.61 that semiconductors and even 
superconductors follow the same scaling law. 

The rest of this section will review the nonlocal response 
functions of the conductor and present calculations that confirm 
these arguments. 



2.3 IVIodel of a nonlocal metal 

In the nonlocal regime, the current-field relationship of the bulk, 
i.e. the conductivity or the dielectric tensor, depends on the 
wave vector and it is necessary to distinguish longitudinal and 
transverse response functions. 
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Figure 2. Reflection coefficients Im Vs.p vs. wave vector in the 
evanescent sector, a) s-polarization and b) p-polarization. The curves 
show specular and diffuse scattering at the boundary, and bulk re- 
sponses given by the local Drude model (solid), and in the nonlocal 
Boltzmann-Mermin model (dotted). Parameters are for gold (cr — 
2.5 X 10''n-^m-\ I = 18 nm, 7 = 6.7 x lO" s"^). Frequency 
uj/2tv = 1 MHz. 



Very basic descriptions taking into account some nonlocal 
effects are hydrodynamic models, see, e.g., Ref. | T3p6| . These 
approaches are valid in a restricted momentum range. The sim- 
plest model does not lead to any change in the s-polarization 
with respect to the local limit, unless some some phenomeno- 
logical transverse response is introduced. A more substantial 
description of a metal is given by the Boltzmann-Mermin (BM) 
model 



cr;(w,k) = - 
0-4(0;, k) 



70- 3uju'^fi{0,u) 
7- iwcj + ijfi{0,u) 

7(7 



-ftiO,u), 



(8) 
(9) 



7 — ICJ 

where a is the DC conductivity in the local limit, 7 describes 
the broadening of the electronic states at the Fermi level due to 
scattering. The variable u contains the momentum dependence 
via the dimensionless functions 

+ 1" 



Mo,u) 



1 
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3 3 



Din 



kvp 
f 1" 



u-1 



(10) 
(11) 



where vp is the Fermi velocity, and the logarithm is taken with 
a branch cut along the negative real axis. These expressions are 



obtained in the limit k/2kp ^ 1 from a more general model 
due to Lindhard fTOl, hence the redundant first argument 0. 
This assumption is reasonable because for our purposes, the 
relevant wave vectors are in the range k ^ 1/z, much smaller 
than the Fermi momentum kp sa 10^^" m~^. 

The relevant frequencies for magnetic transitions, typically 
in the rf- to microwave range, lie in the Hagen-Rubens regime 
w ^ 7 ^ Wp where ujp is the conductor's plasma frequency. 
In this case, spatial dispersion is obviously encoded by the pa- 
rameter u w i/fc£, and the mean free path £ — vp/'j sets the 
relevant scale. In the local limit k£ <^ 1, both conductivities 
reduce to the Drude form 



cri,t(a;,k) -> cr(iLj) = 



Wp£o 



7- 



iw 



(12) 



This is the regime of the normal skin effect where Eq. (|7]) ap- 
plies. 

While the conductivity describes the bulk response of the 
conductor, the specific properties of the surface have an im- 
pact on the nonlocality of response, too. The calculation of the 
reflectivities requires the solution of the electromagnetic scat- 
tering problem at the surface. If the bulk conductivity depends 
on the wave vector, Fresnel's equations do not hold any more, 
and one has to introduce additional boundary conditions for the 
current density at the inner surface. The latter are modeling the 
way charge carriers are scattered there. 

The simplest assumption is that of specular reflection of 
charge carriers [A 7.il0, 13J. Diffuse scattering 1 11 131 or a 
general combination of both mechanisms can be included, but 
must be treated with care to ensure that charge conservation 
holds at the surface |13|. Severe as it may be, this problem 
only occurs when electiic fields have components perpendicu- 
lar to the surface, i.e. in the p-polarization, while the nonlocal 
effects considered in this work involve the s-polarization. In 
addition, it is well known that the scattering mechanisms give 
little differences for the anomalous skin effect [4J. M uch larger 
corrections occur, e.g. due to surface roughness |57[. 

The calculation of reflection coefficients at the surface of 
a nonlocal metal is described in AppendixJB] The resulting re- 
flection amplitudes are shown in Fig. [2] for both polarizations. 
We plot the absorption Im r^p which is proportional, by reci- 
procity, to the radiated noise power. In the s-polarization, all 
models converge to the local scenario as p ^ l/£, as ex- 
pected. Spatial dispersion leads to reduced noise for wave vec- 
tors p ^ l/£. The impact of nonlocality is much more im- 
portant in p-polarization. The increase of p-polarized absorp- 
tion in the range l/£ < p < 1/A has been discussed previ- 
ously |1021,44,58|;itis due to the internal photo-effect (cre- 
ation of particle-hole pairs. Landau damping). The results with 
the diffuse boundary condition introduced in Ref. |fT2) deviate 
from the local limit already in the range p > 1/S, and become 
independent of the bulk conductivity. This is likely to be an ar- 
tifact due to the violation of charge conservation, as discussed 
inRefs. |[n}{l4). 

Note that the factor (w/ck)^ ^ 10^^^ is very small where 
the p-polarized absorption peaks so that it is a good approx- 
imation to neglect this part in the Green's function (J3]l. This 
polarization is essential, on the contrary, for situations sensi- 
tive to surface charges and electric surface fields |13|, such as 
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heat transport |20 |, heating of trapped ions p2) or the electric 
dipole contribution to dispersion forces ||59). 



2.4 Short-distance approximation 

We derive here the asymptotic form of Eq. (J6]l, third Hne. Within 
the approximations introduced above, the Green tensor Q can 
be calculated from 



the normal skin effect was neglected. We conclude that for the 
miniaturization of atom chip experiments, a large mean free 
path is advantageous. Crystalline metals may push i into a 
range that is achievable with atom chip traps. The other pos- 
sibility may be chips based on pure semiconductor substrates 
that we discuss now. 



2.5 Semiconductors 
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The distance range is now z ^ £so that the relevant wave vec- 
tor range is p ^ 1/z ^ 1/i. We start from an expansion of 
the Boltzmann-Mermin model ([8]l-(|9]) at small values of u. The 
limiting form of the reflection coefficients for specular scatter- 
ing of charge carriers is found as (see Appendix IB]) 



p^ 1/i, spec, scatt.: Imrs(p, w) 



1 



p36^{uj)e 



(14) 



while the diffuse boundary condition yields a result smaller by 
a factor 3/4. This is similar to the findings of Reuter and Sond- 
heimer pi for the anomalous skin effect. The following anal- 
ysis assumes specular scattering. The power law of Eq. ([14]) il- 
lustrates the reduction of noise by spatial dispersion (the decay 
with momentum p is faster) and agrees well with a numerical 
calculation, as illustrated in Fig. [2] 

We split the integration range at p ^ l/£ and replace for 
p ^ 1/i the reflection coefficient rs{u!,p) by the nonlocal ap- 
proximation ( 14 1. The integral then gives 



oo 

/ dpp'^rs{uj,p)e'^ 



2pz 



5'^{uj)(. Jl/f. P 



dp _ 



2pz 



Ei{2z/e) \n{£/2z)-jE 



6^{uj)£ 



d^{uj)£ 



(15) 



Here, Ei(z) — J dtt^^e^^* is the exponential integral and 
7b « 0.577 is the Euler-Mascheroni constant. In the range 
p ^ l/£, the reflection coefficient rs{uj,p) is approximately 
equal to its local form [Eqs.(32 33]l]. The integral then gives 
(we assume £ -^ 6) 



The previous analysis can be generalized to other classes of 
conducting materials. Strongly doped semiconductors (where 
the electron gas is degenerate like in a metal) may be described 
by the BM model ((8|)-(fTT|) with modifications only in the values 
of the parameters. For frequencies well below the gap, a back- 
ground dielectric constant et appears due to the static interband 
polarizability, but this does not play a role for magnetic near 
field noise. Weakly doped semiconductors (non-degenerate elec- 
tron gas), on the other hand, are not ruled by Fermi statistics, 
but by a thermal distribution with the characteristic velocity 
vt = {3kBT/my^^ taking over the role of the Fermi velocity 
vp. Here, Eqs. ( [TO] i and ( fTTj i should be replaced by |J60]l6T] 



ft{0,u) 



1 



Imo,u) 



dz^^Pt^ . (18) 



We observe that this is numerically very close to Eqs. ( [TO] i, ( fTTj ) 
in the local regime k£ <^ 1 and differs only by a numerical 
factor in the deeply nonlocal regime. Up to this changed pref- 
actor, the preceding calculations for the Vg coefficient and the 
Green's tensor carry through so that the physics is qualitatively 
the same. 

Let us consider typical numbers that can be found from ex- 
periments on charge transport in silicon ||62J. An n-type semi- 
conductor with a rather low doping of 3 x 10^^ cm~^ is char- 
acterized by a skin depth S > 1 cm in the MHz range. At room 
temperature, vt ~ 1-2 x 10^ m/s and £ w 60 nm, quite com- 
parable to the value for gold. Since it is proportional to the 
DC conductivity, the magnetic near field spectrum is smaller 
by orders of magnitude compared to a metal. Cooling the sam- 
ple down to 100 K reduces phonon excitations and enhances 
the conductivity. Yet, the thermal velocity drops also to vt — 
6.7 X lO'^m/s, so that the mean free path is barely larger, 
^ sa 85 nm. At even lower temperatures the conduction band 
occupation freezes out and the response becomes local. 



p dplmrs{uj,p)e 



-2pz 



l + 0{£/S) 
25^{u)£ 



(16) 



Summing the two contributions gives the approximate Green's 
tensor (always for z <^ £) 



ImHij{z,uj) 



HqIS^j + ZiZj 



In 



£ 
2z 



■IE 



8TrS^{uj)£ 

(17) 
This is the third regime of the Green's tensor (J6]l discussed 
in Sec. 12. 21 and corroborates the statement that a nonlocal de- 
scription predicts less noise at short distances compared to a 
local one. We have thus generalized a similar result reported in 
Ref. p5| within the context of nuclear spin relaxation, where 



2.6 Superconductors 

The nonlocal response of superconductors in the microwave 



range was discussed in Refs. |64 65 1. It was argued by Rick- 
ayzen |66| that screening in a superconductor does not differ 
greatly from a normal metal. This is because all charge carri- 
ers contribute to screening, while the specific properties of a 
superconductor are determined by the states close to the Fermi 
level. We thus expect the noise spectrum to be characterized 
by the same logarithmic asymptote ( [T5| ) found earlier This is 
indeed confirmed by numerical calculations of magnetic noise 
near a niobium surface, the results of which are shown in Fig. 
|3] We have evaluated the non-local BCS conductivity within 
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Figure 3. Spin flip (loss) rate near a surface made of supercon- 
ducting niobium (solid curves) and a fictitious normal metal (dashed 
curves). Nonlocal effects are given for both cases as compared to the 
local limit. Parameters for niobium at T = 0.5Tc follow Ref. |63| 
ujp = 1.33 X 10^''s~\7 = 2.38 x 10""s"\rc = 9.2K, Z\(0) = 
1.9 kBTc ^ h X 7.9 GHz, £ = 5.7 nm. We chose the Larmor fre- 
quency uj/2tv — 1 MHz and 6 — 0. Note that Ref. (64| gives a mean 
free path larger by a factor of 4 and a smaller value of 7. For better 
comparison the normal metal is obtained by closing the superconduct- 
ing energy gap and rescaling the plasma frequency. Gray dashed lines 
indicate the asymptotes to the superconductor. Losses due to the free 
space black body spectrum are not visible on this scale. 



the approach of Mattis and Bardeen including disorder scatter- 
ing | [67) , using the expressions of Popel |64|. The local limit 
recovers correctly the results from Refs. |68 69|. For compar- 
ison we also give the curves for a fictitious normal metal, with 
the same parameters except that the gap is closed {A{0) = 0). 
The description then coincides with the nonlocal BM model. 
The metallic plasma frequency was adjusted in order to take 
into account the redistribution of the spectral weight by disor- 
der, as discussed in Refs. |69,70|. 

In the local regime z ^ i, the superconductor indeed shows 
a strong reduction of magnetic noise. This happens because 
the relevant frequency is below the gap, fku <C ^(r), and 
magnetic fields are well screened by the Meissner effect. At 
large distances, we find good agreement with the expressions 
for the lossrate obtained by Skagerstam et al. [71 1 for a two- 
fluid model, if the same re-scaled plasma frequency is taken 
into account (gray dashed asymptote in Fig.[3]l. In terms of the 
Green's tensor. 



ImHxxiz^uj) 



MoAi(T) 

4:nS^{uj,T)z* 



(19) 



The Meissner-London length Xl{T) = {c/ujp)A{0)/A{T) sa 
{c/(jjp)[l — (r/Tc)^]^^ determines the penetration depth for 
quasistatic fields, while the skin depth d{uj, T) involves the 
conductivity for the normal fluid fraction. The Meissner ef- 
fect becomes inefficient, however, if the spatial scale 1/pof the 
noise field becomes comparable or smaller than the penetration 
depth Al(T). The loss rate for z ^ Xl{T) then approaches 
the 1/z asymptote of a normal conductor, see Fig. [3] At shorter 
length scales z <C ^, we recover the logarithmic scaling law 
found before for the normal conductor (thin dashed line). This 
illustrates the very general character of this regime that does 
not depend greatly on the material class. 



3 Lateral coherence 

A nonlocal conductivity creates spatial correlations in the cur- 
rent fluctuations below the surface which reduce the overall 
magnetic noise level. It is to be expected that this leaves also a 
signature in the correlations of the field. These correlations are 
universal for blackbody radiation |36 1 and have been studied in 
Refs. I38j-40J for homogeneous media and in Refs. 141,-461 for 
the near field of bodies. The spatial correlation length can be 
much larger or much smaller than the wavelength, depending 
on the polariton modes that dominate the electromagnetic field 
noise. We find in this section that the correlation length is con- 
nected to the mean free path as a direct consequence of the bal- 
listic motion of the charge carriers. This should be contrasted 
to electric fields near nonlocal solids where the spatial corre- 
lations were found to differ from the local description only at 
distances comparable to the Thomas-Fermi length A, with the 
mean free path £ not playing any role f44l. 

We are interested in the correlation between fields at a fixed 
height z from the surface and laterally separated by a distance 
p. We define the coherence function as the cross-correlation 
spectrum of the normally ordered field operators in frequency 
space: 



B^j{p,Z,U!) 



dt{:B,{p,z,t)Bj{Q,z,0) 



Aujt 



(20) 



The fluctuation-dissipation theorem for normally ordered op- 
erator products [72J provides the link to the two-point Green's 
tensor 



Bij [p, z, cj) = 2hn Im "Hy (r , r', uj) 



(21) 



where the points r, r' are located at the same height and later- 
ally separated by p. This generalizes Eq. (|2]l. A general integral 
form is given in Appendix[A] We drop the frequency arguments 
for simplicity in the following. 



3.1 Local limit 

Where a local description of the metal is sufficient, the co- 
herence function in the near field can be evaluated asymptot- 
ically by expanding all integrands for values p ^ 1/(5, w/c. 
The resulting integrals have the form /^ dp e~''^P^ Jn{pp)p'"' , 
(n, m e No) and can be evaluated exactly. 

We find that at distances (. <^ z <i 6, the magnetic co- 
herence tensor depends on the distance d = ^J p^ + (2z)-^ be- 
tween one observation point and the mirror image of the other. 
The tensor elements are very well approximated by (see Fig. [4]) 



^y(p,2) 



4^2 S, 




(22) 



Here, the noise spectrum B^x (0, z) for p = was introduced as 
a convenient scale [see Eq. (J6])]. The axes are chosen such that 
the X-axis points along the separation between the two obser- 
vation points. These quantities are independent of the specific 
material properties and depend only on the ratio p/z, i.e. the 
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geometry of the system. An equivalent form for the xx compo- 
nent was already given in Eq. (33) of Ref. f28l, see also Ref. 
||2|. Note that the xz cross-correlation was missed in Ref. |45 1. 
The coherence functions decay on a typical length scale. 
For example, the zz-component (and similarly for the other 
ones) is characterized by the correlation length 



i<t:z 



A 



loc 

zz 



2\/3; 



(23) 



where Bzz{p, z) drops to half its value at p = (see Fig|4k)). 
The approximate forms of Eq. (J22]i are not valid far beyond the 
correlation length where some correlation functions become 
negative, as shown in Fig.pk). The agreement with the asymp- 
totes is so high, however, that the curves are hardly distinguish- 
able. The peak in the crossed xz-correlation arises from light 
paths that are reflected from the surface at oblique angles and 
whose fields are polarized in the xz-plane. 
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Figure 4. Components of the coherence tensor of the magnetic field 
near a gold surface. Solid (dotted) lines: Boltzmann-Mermin model 
with specular boundary condition (Drude model), a) Distance z — 
1 fim much larger than mean free path ^ = 18 nm: all curves coin- 
cide with the local approximation p.2\ (dashed gray curves), b) Non- 
locality becomes important at short distances. For better visibility we 
set ^ = 10 ^m S> z = 10 nm by lowering the scattering rate 7. The 
thick dashed line gives the rather coarse approximation from Eq. j28[ > 
(dashed gray curve). Frequency iu/2n = 1 MHz. 



3.2 Nonlocal metal 

In the near field z <^ £, nonlocality leads to a larger coherence 
length than predicted by the local scenario, as is clearly visible 
in Fig.|4]5). 



hi this section we extract the relevant scales by evaluating 
the two-point Green's tensor Hzz {p, z). It is necessary to inter- 
polate the reflection coefficient rs{p) between the local and the 
nonlocal limits to avoid an unphysical logarithmic divergence 
at the lower bound. The simplest choice is the Pade approxima- 
tion 



Imrs(p) 



1 



I + pi/2 "^ WpHiTplJi) 



(24) 



where the last form is appropriate for p':$> 1/5 and the small -p 
divergence is removed by the factor p'^ under the integral [see. 



e.g. Eq. (13 1]. The integrals then give 
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(25) 
(26) 



where Ei is the exponential integral defined after Eq. ( [TS] ). We 
split the integral atp ~ xo/p where xq « 2.356 is the first zero 
of the Bessel function Jo{x). For < p < Xq/p, Jo{pp) is 
replaced by a spline j{x) = 1 + bx^ + ax^, and replaced by its 
asymptote j{x) ~ Y^2/7ra;cos(7r/4 — x) for larger arguments. 
It turns out that the first interval gives the dominant contribu- 
tion, since the oscillations beyond xo/p provide a cut-off for 
the integrand. We obtain 



'Hzz{p,z) ~ HzziO^z) 
bp 



Pa 



167r(52^ 



-4Ei(2a;oz/p) 



+ ^(l-(l + 22^o^/p)e-2-«^/'') 



ap 



(1 - (1 + 2x,,z/p + 2(x^z/pf) c-2-o./p) 



-°(?4) 



(27) 



A careful glance at this expression shows that in the regime 
z ^ p ^ I , the first line dominates. The decay of the lateral 
coherence is therefore logarithmically (see Fig.|4]i, as the small- 
argument approximation to the exponential integral illustrates 
(thick dashed line in FigpV))) 



B^^{Q,z) 



2E^{2xqz/p) 
Ei(4z/^) 



In [{2xQz/pf 
In [4z/^] 



This yields a coherence length 



(28) 



(29) 



much larger than its local counterpart (23 1. The spatial corre 



lations of magnetic near fields are therefore linked to the char- 
acteristic mean free path of ballistic transport. We recall that 
a similar discussion in Ref. | [44) for electric correlations did 
not find Eq. (29 1 involving the mean free path I, but rather the 
Thomas-Fermi screening length A which is typically smaller. 
The impact of spatial dispersion is therefore somewhat easier 
to reveal by analyzing magnetic fields. 
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4 Discussion 

We have found that the scattering mean free path of charge car- 
riers sets the distance scale for the onset of nonlocal effects in 
the near field noise. For evanescent modes that dominate the 
near field, the fraction of the metallic volume that contributes 
to noise is limited by screening more efficiently than by the 
normal skin effect. The present calculation indicates, therefore, 
that loss rates are actually lower at short distances than pre- 
dicted by local conductivity models (Ohm's law). This noise 
reduction sets in at atom-surface separations comparable to or 
below the mean free path. In a clean (crystalline) metal t may 
take values in the order of 1 /j,m which is at the limits of the ex- 
perimentally accessible region (cf. the gray box in Fig.[T]l: typ- 
ical traps operate at distances of 1 . . . 100 /im from the surface 
and can resolve lifetimes up to lO^'^ ...Is. The data shown in 
Fig. [T] are calculated at room temperature. Since both the con- 
ductivity and the mean free path depend on temperature, it is 
worth investigating whether the lifetimes of magnetic levels 
may be tuned by cooling the atom chip device. This strategy 
is hitting a limit in the extreme near field (distance z <C t): 
the noise power becomes independent of the scattering rate of 
carriers (the Drude parameter 7), and also the details of the 
scattering mechanism of charge carriers at the inner surface 
become irrelevant. The near field spectrum [given in Eq. ( [T5] l] 
has a rather general character and is expected to apply to doped 
semiconductors and even to superconductors, as our numerical 
calculations show (Fig.l3]l. 

We have also analyzed the asymptotic form of the noise 
correlations in the short-distance range and found that mag- 
netic fields are laterally coherent on a scale ^ ^(£ + z)z. This 
implies for an atom chip environment that fluctuating forces 
due to magnetic field gradients are smaller (their spectral den- 
sity scales roughly with the inverse square of the correlation 
length). Also when matter-wave interferometry involves the spa- 
tial splitting of a thermal cloud or condensate, the increase in 
spatial coherence makes the device more robust against deco- 
herence from magnetic noise (see | p8]|73][74| ). The increase 
in spatial coherence intimately relates to the reduction of heat 
transfer via fluctuating near fields because the effective number 
of channels is inversely proportional to the "coherence area". 
For a more detailed discussion of this link, see Refs. | |75][76l . 

Patch potentials due to adsorbates on the surface are known 
to add significantly to the electric field noise, relevant for ion 
traps |77, 781 and systems involving precisely tuned electric 
dipole transitions, such as Rydberg states | [79|j8T[ . Yet they 
will not effect the magnetic case. Static patches have no impact 
on the magnetic noise spectrum and magnetic surface-dipoles 
due to adsorbed atoms have only minor impact: A static charge 
trapped at a distance of 1 /im from an atomic-scale electric 
dipole results in an interaction energy in the order of 20 MHz, 
while the interaction between two such dipoles gives only 1 kHz, 
and the magnetic counterpart for two magnetic dipoles of one 
Bohr magneton /i^ at the same distance gives 1 mHz. More 
prominent sources of magnetic noise might involve diffusive 
currents confined to a surface layer. Their effect may still be 
negligible, however, as the analysis of Ref. |82| has found. 

All of these results imply that nonlocality may be visible 
at the edge of what is feasible with atom chips. Still, operating 
a chip trap at short distances is fundamentally limited by the 



Casimir-Polder interaction that deforms and breaks the trap- 
ping potentials. Alternative setups might, therefore, address the 
broadening of magnetic transitions spectroscopically, e.g. us- 
ing evanescent-wave based surface traps as in Ref |83 1 or op- 
tical tweezers. One may also think of muonic or nuclear mag- 
netic moments, as used in the experiment of Ref. [84) on spa- 
tial dispersion in superconductors. 
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A The magnetic Green's tensor 

The Green's tensor describes the magnetic field 

S,(r,w)=H„(r,r',w)/x,(w) (30) 

radiated by a point-like magnetic dipole source placed in r'. 
Near a single surface, this field consists of a free-space part 
and a reflected contribution, and therefore H = 1^. + H . 

For coinciding spatial arguments, the imaginary part of the 
free space Green's tensor is given by 



'^-"^--'^o^.^. 



lmHi(r,r,a;) 



Gttc'^ 



(31) 



A detailed discussion including regularization procedures is 
given in Refs. 1 85 86). The general expression for the reflected 
Green's tensor reads 

00 



where n ~ \J]P- — iJ^ jc?. Re k > 0, Im k < 0, is the propa- 
gation constant. The tensor structure is included in 

'«'(^o(w)-^) -V^^UVP) 
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where J„(x) indicate Bessel functions of the first kind. The 
relative separation along the surface has length p = [|r — r' p — 
|z — z'p]^'^ and points along the x-axis. In the limit r — > r', we 
obtain the one-point reflected Green's tensor (|3]l that depends 
only on the distance from the surface. It is easy to check that the 
free-space contribution (31 1 is typically negligible at distances 
from a conducting surface smaller than the vacuum wavelength 
c/uj, and we can drop the superscript ^. 
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B Nonlocal reflectivities 

B.1 Reflection coefficients 

For p, s-polarized incident waves, reflectivities are given by 
lITOl 






Zs + ZO 



(32) 



where the surface impedances Zg, Zp are made dimensionless 

by normalizing to the impedance of free space, Z = yT/o/eo- 
The well known results from local theory are ( 19 87) 



^loc ^ ^/s{0J)uJ^/c^ ~ p^ ^lo^ ^ 



uj/c 



e{uj)uj/c 



\J e{LS)ijj^ I (? — ]p- 



(33) 



where e(uj) is the dielectric function (we set liiiS) — 1), and 
the reflectivities ( (32] i reduce to Fresnel formulas. The surface 
impedance on the vacuum side, Z" Z", is obtained by setting 
e(w) = 1. Our focus is on metals at low frequencies where 
in the local limit, e{oj)uP' /c^ ~ 2i/5'^{ijj) in terms of the skin 
depth (J7|, and on the sub-wavelength limit p '3> uj/c. This gives 
an r, coefficient 



ip(5(cj) + v/2i-p2(52(w) 

The large-momentum asymptote (p ^ 1/(5) is rs{uj.p) -^ 
i/{2p^6^). 



B.2 Impedances for specular and diffuse scattering 

The additional boundary condition that charge carriers undergo 
specular reflection at the inner metal surface can be exploited to 
extend the conducting half-space into a fictitious homogeneous 
medium, using a similar symmetry for the fields pll? 10 1. The 
resulting impedances are 



zr'iLo,p) 



zr%c^,p) 



2iLj 
en 
2iuj 



dq 
dq 

fc2 



et{uj, k)w2/c2 — fc2 



e/(iu,k)w2/c2 



(35) 



(36) 



St{tL),'k)uj^/c^ — fc2 



2 The dielectric 



where flie medium wave vector is fc^ = p^ + q 
functions are et,i = £b + i(Tt,i/£o^ with a background polar- 
ization £{, that drops out in our regime. If the medium is local, 
Et^i do not depend on k, and a direct calculation of the integrals 
brings us back to the Fresnel impedances ( (33] l. 

The diffuse reflection of conduction electrons at the sur- 
face has been considered in Refs. |4,111288|. We have used 
the s-polarized impedance obtained in Ref. [Jlj where the ad- 
ditional boundary condition is implemented via the so-called 
dielectric approximation. This means that in the basic linear 



current-field relation [Eq. ([TJ], the volume integral is restricted 
to the medium-filled half-space alone. Solving a Wiener-Hopf 
equation that follows from the Maxwell equations, one gets an 
impedance 



zr{^,p) 



ITTUJ 
C 



dqlog 



k^ — et{uj,'k)uj^/c^ 



-1 
(37) 



This recovers correctly the local limit at all values of p. How- 
ever, this is not true for the p-polarized impedance given in the 
same work and used in the numerical evaluation in Fig. [2] The 
violation of charge conservation at the surface in the dielectric 
approximation has been discussed in Refs. 1 12 14). We con- 
tinue here with the specular boundary condition. 



B.3 Limiting behavior of the reflectivities 

For a metal with specular reflection of charge carriers, Eq. ( (35| ) 
is rewritten in a dimensionless form by rescaling the wavevec- 
j2 — pt. Making the low-frequency approxi- 



tor k = \fp^ 

mation described before Eq. ( 34 1, the non-local dielectric func- 
tion becomes [see Eqs.(|7]|9]l] 



et(w,k) 



2i/t(0,i/fc^) 



(38) 



For simplicity, we drop in the following the frequency argu- 
ments and the redundant argument of the Lindhard function ft- 
We are left with the integral 
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2uj/c 



(39) 



We expanded the denominator in the sub-skin depth regime 
(p ^ 1/6, l/£) and replaced, in the last step, the transverse 
Lindhard function [Eq. ( [TT] l] by its asymptote for large pt. 

Mi/pit) = — +0[ipit)-^]. (40) 

This result complies with a similar calculation carried out in 
Ref. |_89^|. 

The first term in Eq.([39| corresponds to the free-space sur- 
face impedance in the sub-wavelength limit. The reflection am- 
plitude (|32jl therefore becomes 



spec 



b) 



ZIP^^{p)+iL0/{cp) 



^spcc 



(p) — iuj/{cp) p^S^i 



(41) 



which gives Eq. ( 14 1 



In the scenario where charge carriers are reflected diffusely 
rather than specularly we start from Eq. ( (37| ). The same substi- 
tution of the integration variable gives 



7diff 



(P) 



ITTUJ /(cp) 
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log 
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(42) 
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The integrand is expanded in l/(p(5)^, and the integration can 
be performed explicitly 



dtt 



dtt 
SttI 



log 
log 



i' 



1 



t^-2i{pSt)-^Mi/p£t) 
2iMi/p£t) \ 



2pHH 
The resulting impedance reads 



t2 

if pl > 1 



(pStr J 



7diff 



(P)«- 



iu! 3io/c 

^ ^ 2p^s^e 



(43) 



(44) 



This differs from Eq. ( (39| ) only by a factor 3/4 in the second 
term, so that the reflection amplitude is smaller by this number 
compared to the specular boundary condition, Eq. (pT). 
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